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NUMERICAL  RESULTS  FOR  LOW  FREQUENCY  SCATTERING 
BY  ELLIPTIC  CYLINDERS  AND  BY  ISOLATED 
SEMI-ELLIPTIC  PROTUBERANCES 


by 

J.  E.  Burke,  E.  J.  Christensen,  and  S.  B.  Lyttle 

ABSTRACT 

The  low  frequency  approximations  ("closed  form"  and  series) 
derived  previously  for  the  fields  scattered  by  elliptic  cylinders,  and 
by  semi-elliptic  protuberances,  are  applied  numerically.  For  the  two 
cases,  Ej  or  H  parallel  to  the  generators,  the  results  presented  in¬ 
clude  total  scattering  cross  sections,  forward  and  back  scattered  inten¬ 
sity  and  phase  curves,  and  far-field  scattering  patterns;  for  various 
angles  of  incidence,  various  eccentricities,  and  for  various  values  of 
fcd-  ^  1.1  (where  k  -zrr/S)  and  ZCL  is  the  major  axis  of  the  scatterers). 
Attention  is  restricted  to  the  low  frequency  range  not  covered  by  pub¬ 
lished  tables  of  Mathieu  functions. 


1.  INTRODUCTION 

The  problem  of  scattering  by  a  perfectly  conducting  elliptic  cylin¬ 
der  is  separable  in  elliptic  coordinates  and  the  solution  can  be  represented 

12  3 

as  an  infinite  series  of  periodic  and  radial  Mathieu  functions  9  1  .  How¬ 

ever,  since  Sieger1  original  derivation  of  this  solution  in  1908,  the  only 

4 

detailed  calculations  based  on  it  appear  to  be  those  of  Morse  and  Rubenstein 
for  the  strip  (the  elliptic  cylinder  with  eccentricity  equal  to  1).  In  this  re¬ 
port  the  exact  series  are  used  to  obtain  numerical  values  for  the  scattering 
amplitudes  for  a  family  of  elliptic  cylinders  (ranging  from  strips  to  circles), 
and  for  the  corresponding  semi-elliptic  protuberances  on  ground  planes. 

For  limited  values  of  kd  ^1.0>  where  and  ZA,  is  the  major 

5,  6 

axis  of  the  cylinder,  tables  and  graphs  of  the  Mathieu  functions  can  be 
used  to  evaluate  the  exact  series  for  the  scattering  amplitudes.  In  the  pre¬ 
sent  report  we  consider  the  essentially  complementary  range  fat  ^  1.1  and 
evaluate  the  series  in  "closed  form".  The  closed  forms  are  obtained  by 


initially  truncating  the  exact  series  and  then  using  values  for  the  Mathieu 

7 

functions  obtained  from  known  low  frequency  approximations  .  The  valid¬ 
ity  of  these  forms  is  investigated  by  comparing  them  numerically  with 
"exact"  results  based  on  tabulated  functions.  In  addition  to  considering 
closed  forms,  we  demonstrate  the  utility  of  elementary  series  approxima- 
tions  (in  powers  of  k  )  for  the  scattering  coefficients  and  for  the  scatter¬ 
ing  amplitudes. 

Although  we  work  with  relatively  simple  analytical  expressions, 

their  numerical  application  is  complicated  by  the  fact  that  four  variables 

are  involved;  the  eccentricity,  the  angle  of  incidence,  and  the  angle 

of  observation.  To  make  the  computations  tractable  they  were  performed 

on  a  Burroughs  220  electronic  computer.  The  extensive  numerical  results 

thus  obtained  are  preserved  in  permanent  tables  and  in  punched  cards; 

the  latter  provide  data  in  a  usable  form  for  electronic  computing  programs 

8  9 

on  related  multiple  scattering  problems  .  It  is  not  possible  to  present 
the  complete  set  of  tables  here  and  we  merely  illustrate  their  contents 
through  a  series  of  graphs.  For  cylinders,  and  for  protuberances,  these 
graphs  include  forward  and  back  scattered  intensity  and  phase  curves, 
total  scattering  cross  sections,  and  far-field  scattering  patterns;  for  var¬ 
ious  angles  of  incidence,  various  eccentricities,  and  for  various  values  of 
*  1.1.  A  variety  of  results  are  given  so  as  to  demonstrate  the  overall  de¬ 
pendence  of  the  amplitudes  on  the  four  variables. 

In  the  following,  we  begin  with  a  brief  review  of  the  scattering 
problem  and  the  derivation  of  the  series  solution  by  separating  variables 
in  elliptic  coordinates.  Then  approximations  for  the  Mathieu  functions 
are  introduced,  and  the  accuracy  of  the  resulting  expressions  (closed  form 
and  series)  for  the  amplitudes  is  discussed.  In  the  final  section  the  tables 
that  have  been  compiled  are  described  and  graphs  of  results  taken  from 
them  are  given. 


Z 
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REVIEW  OF  SCATTERING  BY  AN  ELLIPTIC  CYLINDER 


In  two  dimensions  the  scattering  of  a  plane  wave  by  a  cylinder  par¬ 
allel  to  the  axis  is  specified  in  the  region  external  to  the  scatter er  by 
a  solution  of 

(1)  ^  ~  °>  V  ~o>x+c jy,  k~zn/}\, 

satisfying  prescribed  conditions  at  the  surface  of  the  cylinder.  The  solu¬ 
tion  has  the  form 


(2) 

where 

(3) 


no  =  %(o  +•  u(o> 


X(rj  -  e:kTCOA(^-^ 


represents  the  incident  plane  wave,  and  where  the  associated  scattered 
wave  fulfills 

(4)  U(r)  ^  G  ^  V Z/(tnkr )  >  t 

The  "scattering  amplitude”  indicates  the  "far -field"  response 

in  the  direction  *P  to  plane  wave  excitation  of  direction 

O  * 


For  the  elliptic  cylinder  the  boundary  conditions  are  applied  on 
the  surface 


(5) 


y-£  u* 


where  ZCL  and  zb  are  the  major  and  minor  axes  respectively;  see  Fig.  1. 
In  particular,  for  an  incident  electromagnetic  field  specified  by  either 
4,**  or  bY  (the  two  principal  polarizations),  the 

component  of  the  corresponding  total  field  on  the  surface  of  a  perfect 
conductor  satisfies  either 


(6) 

II 

£ 

or 

(7) 

//  Ea  = 

X/ 

**  tL-tr  ■ 
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To  solve  the  problem  posed  for  the  elliptic  cylinder,  one  introduces 
elliptic  coordinates  (%yQ)  defined  by 

(8)  +  =  r&l'f  ^  jrao6h(Z+le)> 

as  illustrated  in  Fig.  2,  the  coordinate  curves  £  =  constant  and  O  =  con¬ 
stant  are  ellipses  and  hyperbolas  respectively,  and  each  curve  has  foci  on 
the  X  axis  at  Then  Eq.  (1)  becomes 


(9) 


+  -^2+- k£(cc*h2z-  co6£&)  r(K,m  =  o  -  k*(a?-b*) , 


and  the  cylinder  (5)  is  specified  by  the  complete  coordinate  surface  %  —  £ 
where 

(10)  0-  —  £  j  b  —  £  fcunh  ^  . 

By  separating  variables,  the  solution  of  the  boundary  value  problem  can 
be  written  as  a  series  of  particular  integrals  of  (9)  of  the  form  ? 

thus 

(id  =  Je 

where  the  constants  oCn  are  determined  by  applying  the  boundary  conditions 

at  r  = 

In  (11)  the  Sn  are  the  periodic  Mathieu  functions  (even  or  odd,  of 
period  ft  or  £,/T  )  and  the  are  the  radial  Mathieu  functions  of  the  first 
and  second  kind  (analogous  to  Bessel  functions  of  the  first  and  second  kind). 
These  functions  are  discussed  in  detail  in  the  literature^*  ^  and  in  the 
following  we  only  state  those  forms  and  properties  needed  to  make  (11) 
explicit. 


For  a  given  value  of  A  t  the  periodic  Mathieu  functions,  which  we 
denote  by  §  (0)  (instead  of  by  $  (h^CC&B)  as  in  reference  11),  form  a 

o  & 

complete  set  of  orthogonal  functions: 


(12) 


r  zn 

\  5,  6 e)So  ce)d.e  =  ot 

J  cm  °n  ? 


l  sJje>^(e)de- 


m  to 

m -n 


? 


J=e>°  • 
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where  the  present  equal  those  of  reference  11  divided  by  2rr.  Their 

Fourier  expansions  are  of  the  form 


(13) 


nB )  AnCfi'),bUin& > 


A*, 


where  2ii  means  that  m  and  n  have  the  same  parity  and  where  the  val¬ 
ues  of  the  coefficients  depend  on  the  normalization  of  the  S •  C&)'  as  in  refer- 

Jm  * 


ence  1 1  we  use 


SaJ@)  -> 


/torn  mB> 


Sure) 

771  yyl 


for  h  ->  O 


(14) 


4/w -g6>>  - <  ■ 


n 


The  present  functions  are  related  to  the  C-  and  S<*  of  references  Z  and 
1  0  through 


(15) 


% • 


'/* 


So  ~  S6T  6?  M0  ) 

cm  n  ,5h/ 


The  radial  functions  of  the  first  kind  )f  may  be  written  as 

series  of  Bessel  functions: 


(16) 


^ S 5) 

m  77- 1 


The  functions  of  the  second  kind,  Ne  follow  from  (16)  by  replacing  the 
Bessel  functions  by  the  corresponding  Neumann  functions.  However,  for 
small  values  of  fl}  the  expansions  involving  products  of  Bessel  and  Neumann 
functions  are  more  convenient: 


(IV) 


^  -  rk>  (*t)  *  (%-) , 


For  large  values  of  £T;  the  and  behave  asymptotically  like 

their  Bessel  function  counterparts,  i.  e.  ,  JT  (r) J  (fer\  A j  ^  a /  (f?r  ) 

grn  rn  *0  gyn?'  /• 

In  particular,  the  linear  combinations 


(18) 
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are  analogous  to  the  Hankel  functions  and  correspond  to  outgoing  cylindri¬ 
cal  waves  as 


(19)  H§JV -L-nJ\z/(cnhoc*h%)  ~ir**te/Cinkr)  &ikr,  00 • 

In  order  to  express  the  formal  representation  (11)  for  V7  in  terms 

of  the  above  functions,  it  is  convenient  to  consider  the  wave  functions  3^ 

and  il  separately.  Thus,  employing  the  contracted  notation  F 

2  22  Jn  n-oen  ) 

the  usual  expansion  *  of  the  plane  wave  may  be  written  as 


(20) 


m>e)  -  elkrcoA(<p-V  (e)j. COM:1 . 

?  Jn  Jyi  jyl  Jn 


Similarly,  for  both  polarizations,  the  representation  of  u  in  terms  of 
outgoing  radial  functions  is  of  the  form 


(21) 


ufce)  =£jir>a.SX‘%)  5. 

^  Jn  Jn  Jn  Jn  Jyi 


The  scattering  coefficients  d:  determined  by  applying  the  boundary  con- 

Jyl  ' 

dition  (6)  or  (7)  are  given  by 


(22) 

_  m  A- 

*  H-CZ)  *  3 

Jn 

and 

(23) 

*  Jn3 

v£;. 


j. 

Jn 


dX'. 


3? 


^  I  f  e±a.  C&r  H=  . 


The  Mathieu  function  series  representation  for  the  scattering  amp¬ 
litude  follows  by  letting  %  or  kv  become  infinite  in  u  of  (21).  Thus, 
using  the  asymptotic  form  of  tie  of  (19),  we  obtain 


(24) 
where 

(25) 


r- 


^Z/ (inhcakfi  £  ) ^(B/Pq  )  %  ) 


The  amplitudes  for  the  elliptic  cylinders  can  be  used  to  obtain  the 
analogous  ones  for  semi-elliptic  protuberances  on  ground  planes;  see  Fig.  3. 
These  follow  by  taking  twice  the  symmetric  or  anti -symmetric  components 
of  the  cylinder  results  with  respect  to  reflection  of  one  angle  in  the  plane  of 
one  axis;  a  procedure  equivalent  to  that  used  originally  by  Rayleigh^  to 
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(b) 


Fig.  3.  Geometry  for  scattering  by  a  semi-elliptic  protuberance  on  a 

ground  plane.  Case  (a)  shows  the  major  axis  perpendicular  to  the 
ground  plane  (6  =  T]/^  =  h/a  =  p  <  1)  and  case  (b)  shows  the 
minor  axis  perpendicular  (6  =  T|/ 5  =  a/b  =  p-1  >  1);  6  =  0 
corresponds  to  perpendicular  strips,  6  =  1  to  semi-circles, 
and  6  =  oo  to  flat  strips. 
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treat  the  semi-circular  protuberance.  Thus,  if  the  minor  axis  is  in  the 
ground  plane  X  —  O  (see  Fig.  3a),  if  7t-  g )  is  the  direction 

of  incidence,  and  if  H  is  parallel  to  the  generators,  we  obtain 

(26) 

Similarly  for  JE  parallel 

(27)  n-40  %) . 

The  corresponding  results  for  the  major  axis  in  the  ground  plane  (see  Fig.  3b) 
follow  from  (26)  and  (27)  by  replacing  )  by  ^  ?  or  > 

equivalently,  by  interchanging  Cl  and  b , 


The  real  part  of  the  scattering  amplitude  for  the  special  value  ^ 
is  proportional  to  the  total  scattering  cross  section  £).Thus,  for  cylinders 

(28)  Q±  “ 

and  for  protuberances 

(29)  gu  -  -  f ' «=  i-A*-®  -  ■ 

~rr/£- 

For  limited  values  of  ka  >  1  numerical  values  for  the  scattering 

5,6 

amplitudes  (25)  to  (27)  can  be  obtained  by  using  tables  and  graphs  *  of 

the  Mathieu  functions.  If  k&L  is  large  enoughqone  can  use  asymptotic 

13,  14 

methods  to  obtain  high  frequency  approximations  .  On  the  other  hand, 

for  k  <  1.0  the  required  values  of  the  Mathieu  functions  are  not  readily 
available.  They  could  be  calculated  by  employing  tables^'  of  the  Fourier 
coefficients  Sm.  However,  these  are  given  for  values  of  k  -  (&&)(!- 

and  consequently  kcu  and  the  ratio  p  =  b /&  cannot  be  varied  independently. 
In  addition,  this  coupling  complicates  the  calculation  of  the  radial  functions; 
in  general,  for  the  given  h  the  values  of  the  Bessel  and  Neumann  functions 
needed  in  (16),  (17),  etc.  ,  are  not  given  in  published  tables. 


In  the  following  sections  we  restrict  attention  to  computations  in 
the  low  frequency  range  not  covered  by  published  tables.  The  coupling 
described  above  is  avoided  by  using  truncated  series  approximations  for 
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the  Mathieu  functions.  Thus,  in  terms  of  f  =  b/&}  %  —  } and  L—/n  - 

7 

(where  /n  —  0,5772/5“  *  is  Euler's  constant),  for  the  radial  functions  we  use 

*.  =  *£'xi  =  J/^b)  ^  =  £(w-Pv  g«-6?a+?*)] 


(30) 


mc^*° = ic^b)  j = ^  r 

3  <°  J>  e™  (m-l)t  ^ic  ™  L  *  J  * 

pX 


f  &■  c. 


^  >  [M/£] = m/£-  > 

Ne~  4&L+XV  jzv+srw^se3*^*) <•  % 6- $■?-?*+■  jtf-e¥-fer+fjl, 

At,*  xjjfejjf*  fOr^-8(l+?)L)i-  *^3^eV^ 

*31^  "  ■fc01+&s.?+3^*i-2qs-h?+^  =■  NCi(a-<^bl 

^ ^ [1+  *  6 - n , 

where  & +~>  b  indicates  the  interchange  of  ^  and 

7 

Similarly,  for  the  derivatives  with  respect  to  ^  we  use 

,J*' "  *  iV/opV??*;^  -5'-  xpf- 

<3\  -  *V£- £c«/A  < -  |W>  -  £&««,  £-  £fV-  £e<«p  V>/, 


and 


/v„' -  !{/+  ^6+?*~*?L)-r  § ih-+z<?*4-<?4+^p(i+?*)Ql} 

(33)  <  =  ££3+/*p-y*-8^/+e>i- 1- 

+  J£(7tffC$+lW<?S+6& 

N‘ H-  ^  6^,  . 
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Many  truncated  series  approximations  for  the  periodic  Mathieu 
functions  appear  in  the  literature  *  .  For  present  purposes  it  is 


convenient  to  use  those  given  in  reference  10: 
,  'M  1-  4  /  —  l  ^ 


"2  -  "  \4/ 


(34) 


^H£)  ^(m^-  ml 

-(A)  A>Ca&'7<P  <206  S cp  <206  5<P  ,  Cc6<P  ? 

\Z/  l  q.£L0  <?.£!'  3-2 *<>  £9  j  p 

V^V1-  fim-  £,****} , 

\un.y~-  A'M‘r~ Kip" 3f Y - *§£$ 

../AY  (AM.?*?  .  6Mt6±  AirtStP  6pn<P  ) 

Uy  2  ?.Exo  ?.£?  5£J-a~  £9  J, 

z. 


Zj  l 


£8%  5 ; 


/6M  <f  ? 


W -^-(1)  Hr*-  ~ rs 

These  give  most  directly  the  angular  factors  which 

appear  in  (25).  Using  Eqs.  (30)  to  (34)  the  exact  series  for  the  amplitudes 
can  be  evaluated  in  closed  form. 

Explicit  series  approximations  for  the  scattering  amplitudes  have 
been  obtained7  by  inserting  (30)  to  (34)  into  (25)  and  expanding  in  powers 
of  k  (treating  L  as  a  constant).  The  amplitudes  correct  to  order 
thus  obtained  are  given  in  reference  7  and  they  will  not  be  repeated  here. 
However,  for  future  reference  we  note  some  leading  approximations.  Thus 
for  cylinders  with  H  parallel  the  real  and  imaginary  parts  to  order  i^are 

Re-4+(eP,cP')=-tt 

P-  8?L) 


(35) 


V 

-h  £pCoQ64)QQ6 5*%  +cofi<f0 Qo,6 5 A3) 


-  (H£}1 

•f 


s 

S^-t-Aon^A# >25<a]- (7-2p  +  '/?L)&#>‘:PcG6‘g 

j£ 


(7% -£ +4OA0at(’&M‘‘fl  +  fs P004 2<faas£tpo + (L+fi/bm ZV/t'px £%J j 
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Similarly,  letting  P PL^)  ,  the  amplitude  for  £.  parallel  is  speci- 
fied  to  order  fe  by 

=-o-h  §(&j[£a-D)+a-f)fa +cc#*p)j , 

(36)  ^f-m)  -  ^p- -  (jr/[--Dx~£D)  +  ^a-?&Xcc*z<P+  coazcpb  ) 

+n(i+%)Coc&‘p(&^oi-p&wz^6^‘Ji J  . 

3.  ACCURACY  OF  THE  APPROXIMATIONS 


In  this  section  we  investigate  the  numerical  validity  of  the  closed 
forms  and  series  approximations  for  treating  elliptic  cylinders  with  arbit¬ 
rary  eccentricity  and  kd<  1.0 .  The  scattering  coefficients  are  considered 
first  (for  kcL<l.Z~  )  and  then  the  scattering  amplitudes  (for  fed  1.1  ).  In 
each  case  the  range  of  kd  includes  some  values  in  common  with  the  tables 
of  references  5  and  6. 


Consider  initially  the  closed  form  approximations  ^S/72/Z^  -  Cs>^X^5)\C3l)'CCd3fl 


-t 


etc.  ,  for  the  coefficients  &  ' 

on 


for  Yi  =  O  and  1  (the  leading  coefficients). 

For  the  limiting  cases  of  strips  (f  -  o)  and  circles  ='  1^  the  accuracy 

of  these  formulas  follows  by  comparing  them  numerically  with  "exact"  re- 

1 8 

suits;  results  obtained  by  using  tabulated  values  for  the  Bessel  functions, 

14  15 

and  for  the  Fourier  coefficients  *  ,  in  the  series  representations  for  the 

radial  functions  (e.g.  ,  in  Eq.  (16)).  Thus,  as  illustrated  in  Table  I,  for 
p  O  or  L  and  fed  ^ l.£ y  it  is  found  that  the  closed  forms  for  the  leading 
coefficients  practically  equal  exact  values.  Similarly,  for  intermediate 
values  of  p  (i.  e.  ,  0<(p  <  X  ),  and  l<fed<I.£^t he  closed  forms  agree  with 
calculations  based  on  references  5  and  6.  Such  numerical  comparisons  in¬ 
dicate  that  for  the  closed  forms  for  the  scattering  coefficients  are 

uniformly  accurate  for  all  p  and  that  the  approximation  for  Cl£  when  p  =  L 
is  the  least  accurate. 


Numerical  results  for  CL;  with  Yls>xL  show  that  for  the  range  fed<l. 

J/z  y 

of  primary  interest,  and  for  the  larger  range  fed  <  1.1  >  all  but  dy^  and  'Jbn  dj^ 
are  negligible  compared  to  the  leading  coefficients.  These  higher  order 
terms  are  most  significant  for  the  circle,  while  for  the  strip  only  \QmCL' 
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Table  I.  Comparison  of  !,exactn  (E)  and  closed  form  (C)  values  for  the 

imaginary  parts  of  the  leading  non- zero  scattering  coefficients  for 

strips  (p  =  0)  and  for  circles  (p  =  1);  note  that  =  &on  5  0 

for  p  =  0  and  agn  =  a*n  for  p  =  1.0.  The  exact  values  are 

based  on  tabulated  functions  and  the  closed  forms  on  the  formulas 

Im  a7  =  (30)(31)[(30)8  +  (31)8]-1  and  Im  at  =  (32)(33)  [(32)2  +  (33)a]  _1  . 
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l 


needs  to  be  retained.  For  the  circle,  the  sum  of  the  absolute  errors  in¬ 
troduced  by  the  closed  forms  for  is  less  than  2  percent  of  the 

corresponding  component  of  a  leading  coefficient. 

As  shown  in  Fig.  4,  for  certain  ranges  of  the  parameters,  the 
closed  forms  for  the  coefficients  may  be  replaced  by  truncated  series. 

In  general,  the  smaller  the  value  of  p  the  larger  the  range  of  kd  for 
which  the  series  apply;  this  is  illustrated  in  Fig.  4a  by  the  series  for  the 
imaginary  part  of  In  particular,  the  agreement  of  the  different  approx¬ 

imations  in  Fig.  4a,  when  p  O  >  is  typical  of  the  results  for  the  real  and 
imaginary  parts  of  ^€/l  anc*  'in  general,  the  results  for  the  other  co¬ 
efficients  are  not  as  good.  The  maximum  range  of  for  all  p  follows 

by  comparing  the  series  with  exact  results  for  the  circular  cylinder.  Thus, 
Fig.  4b  and  the  p  =  1  curves  of  Fig.  4a  determine  the  maximum  ranges 
of  the  series  for  the  real  and  imaginary  parts  of  the  leading  coefficients 
(  CL~  excepted).  (The  series  for  dZ  are  more  accurate  than  those  con- 
sidered,  e.  g.  ,  the  series  is  valid  for  all  p  if  1.1.) 

The  accuracy  of  (34)  for  calculating  the  angular  factors  SnCcf))SnC%)Mr^ 
follows  by  comparing  them  numerically  with  results  based  on  the  tabulated 
Fourier  coefficients.  Thus,  for  (i.  e.  ,  for  and  l) 

it  is  found  that  the  two  results  differ  by  less  than  1  percent.  These  approxi¬ 
mations  are  least  accurate  for  the  strip  ■=.  o)  for  which  case  h  =■  fe/z/xl  j 
for  the  circle  (p—  Z  O  and  (34)  reduces  to  simple  trigonometric 

functions  (independent  of  fl  ). 

Combining  (34)  with  the  closed  forms  for  the  scattering  coefficients 
leads  to  relatively  accurate  expressions  for  the  leading  terms  of  the  series 
(Z5).  Of  course,  because  of  the  angular  dependence,  the  errors  in  the  indi¬ 
vidual  terms  can  cancel  or  add  up  to  give  more  or  less  accurate  values  for 
the  amplitudes.  As  an  illustration  of  the  overall  accuracy  of  the  final  forms, 
Fig.  5  compares  closed  form  and  exact  results  for  the  circular  cylinder 
for  H  parallel;  thus,  for  example,  the  closed  forms  for  |^.l^are  seen  to 
be  accurate  to  within  3.  5  percent  for  all  angles  of  observation.  Better 
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0  0.2  0,4  0.6  0.8  1.0  1.2 

ka 


Fig.  4(a). 

Fig.  4.  Closed  form  and  series  approximations  versus  ka  <  1.2  for 
(a)  Im  a*  (for  different  values  of  p  =  b/a),  and  for  (b)  the 
real  and  1  imaginary  parts  of  leading  scattering  coefficients  for 
circular  cylinders  (p  =  1).  The  solid  curves  are  closed  forms 

and  the  dashed  ones  are  series  approximations;  the  numbers  give 
the  highest  power  of  k  retained  in  the  series. 
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Fig.  4(b). 
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overall  agreement  is  obtained  for  smaller  f&L  the  circular  cylinder 
with  E  parallel,  and  for  strips  for  both  polarizations;  e,  g.  ,  for  the  case 
of  Fig.  5,  the  closed  forms  and  exact  values  are  practically  identical 
for 

For  special  ranges  of  the  parameters,  the  closed  forms  can  be 
replaced  by  elementary  series  approximations;  e.  g.  ,  (35)  and  (36).  This 
is  illustrated  in  Figs.  6  and  7  which  compare  closed  form  and  series  re¬ 
sults  for  cylinders.  Thus,  for  example,  Fig.  6  shows  that  for  H  paral¬ 
lel,  the  /^^approximation  (35)  applies  (more  or  less)  for  when 

p £0.5-  Similarly,  for  E  parallel,  Fig.  7  demonstrates  that  the  fcf 
series  (36)  can  be  used  for  kd^O.7  when  ^£0.5  The  full  range  of  p 
can  be  covered  by  suitably  restricting  &CL  and/or  the  angles;  e.  g.  ,  the  6* 
series  for  fz  parallel  can  be  used  for  all  p  if  fad  0.7  and  the  direc¬ 
tion  of  observation  is  not  near  the  back  direction  (see  Fig.  7b). 

4.  TABLES  AND  GRAPHS 

The  low  frequency  approximations  described  in  the  previous  sec¬ 
tions  have  been  applied  numerically  in  detail.  Most  of  the  calculations 
were  performed  on  a  Burroughs  220  electronic  computer  and  the  results 
are  preserved  in  permanent  tables.  For  cylinders  and  for  protuberances, 
for  £  and  H  parallel,  these  tables  include  closed  forms  for  the  real  and 

Aj  /V 

imaginary  parts  of  the  scattering  amplitudes,  and  for  the  corresponding 
intensities,  for  the  ranges 

i.i)  A(hL)=o.£j 

(39)  =  b/cLH  l.o9  =&.£5- 

%  -o  ,10°  4-5°  SO°,  90 “ 

5SO°p  =  10°  . 

In  addition,  tables  of  forward  and  back  scattered  intensities  and  phases, 

and  of  been  compiled  for 

^  9£?(Z  j5°<AcP<  10°)  and  the  values  of  p  and  kd  of  (37);  the 
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Fig.  6.  Plots  of  |  g+(cp>  45°)  | 3  versus  cp  for  (a)  fixed  p  =  b/a  =  0.5 
(and  different  values  of  ka),  and  for  (b)  fixed  ka  =  0.  5  (and 
different  values  of  p);  the  dashed  curves  are  based  on  the  k4 
series  (35). 
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Fig.  7.  Plots  of  |g  (<p,  45°)  | 2  versus  cp  for  (a)  fixed  p  =  0,  5  (and 

different  values  of  ka)  and  for  (b)  fixed  ka  =  0.7  (and  different 
values  of  p) ;  the  dashed  curves  are  based  on  the  ka  series  (36). 
For  case  (a)  the  series  equals  the  closed  form  when  ka  <  0.  3. 
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smaller  increments  of  were  used  for  protuberances  for  near -grazing 

incidence  (with  respect  to  the  ground  plane)  to  facilitate  obtaining  the 

9 

analogous  results  for  elliptically  striated  surfaces  . 

It  is  not  possible  to  present  the  complete  set  of  tables  here,  and 
in  the  following  we  merely  illustrate  their  contents  through  a  series  of 
graphs.  These  are  given  in  four  sets  corresponding  to  cylinders  with  £ 
perpendicular  or  parallel  (Set  I  and  Set  II),  and  to  protuberances  with  E 
perpendicular  or  parallel  (Set  III  and  Set  IV).  The  individual  graphs  are 
labeled  Graph  1-1,  Graph  1-2,  etc.  ,  where  the  integer  always  refers  to 
the  same  quantity.  Thus,  without  regard  to  a  specific  set,  the  graphs  are 
described  in  general  by  the  following: 

Graph  1:  k/4  times  the  total  scattering  cross  section  (i.  e.  ,  the  negative 

of  the  real  part  of  the  scattering  amplitude  in  the  forward  direc¬ 
tion)  versus  the  angle  of  incidence  (cp  ),  for  different  values 
of  p  =  b/a  and  ka. 

Graph  2:  Forward  scattered  intensity  versus  the  angle  of  incidence, 
for  different  values  of  p  and  ka. 

Graph  3:  Forward  scattered  intensity  versus  the  angle  of  incidence, 
for  different  values  of  p  and  cpQ. 

Graph  4:  Back  scattered  intensity  versus  the  angle  of  incidence,  for 
different  values  of  p  and.  ka. 

Graph  5:  Back  scattered  intensity  and  phase  curves  versus  ka,  for 
different  values  of  p  and  cpQ. 

Graph  6  and  higher: 

Far-field  scattering  patterns  versus  the  angle  of  observation 
(cp),  for  different  values  of  p  ,  ka,  and  cpQ. 

More  detailed  descriptions  of  the  individual  graphs  (e.  g.  ,  the 
explicit  values  of  the  parameters  considered)  are  listed  at  the  beginning 
of  each  set. 

The  authors  wish  to  thank  Miss  Mary  Brockett  who  drew  the  graphs. 
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Parallel. 


SET  I:  Graphical  Results  for  Elliptic  Cylinders  and 

Graph  I- 1:  -Re  g+(cpQ,  <P0)  versus  90  for  ka  =  5'  9, 

and  1.1,  and  different  values  of  p.  * 

Graph  1-2:  |  g^(cp0>  cpQ)  | 2  versus  cpQ  for  ka  =  0.5,  0.7,  0.9,  and  1.  1 

and  different  values  of  p.  * 

Graph  1-3;  [gt(cp  ,  cp  )|2  and  the  phase  of  g  ,  (cp  ,  m  )  versus  kacl.l, 

+  o  o  *r  o  To  “ 

for  cp  =  0°(a),  45°(b),  and  90°(c),  and  different  values  of  p.  * 

Graph  1-4:  |  g_j_(Tr+cpo>  cpQ)  j  a  versus  cpQ  for  ka  =  0.5,  0.7,  0.9,  and  1.1,  and 

different  values  of  p.  * 

|  g^(ir+cpo,  cpQ)  |  s  and  the  phase  of  g^(‘n’+cp0,  cpQ)  versus  ka<l.l, 
for  =  0°(a),  45°(b),  and  90°(c),  and  different  values  of  p.  * 

| g+(cp,  cpQ)  | *  versus  cp  for  cpQ  =  0(a),  45°(b),  and  90°(c), 
ka  =  0.  3,  and  different  values  of  p.  * 

Graph  1-7:  |  g+(cp,  cpQ)  | 2  versus  cp  for  cpQ  =  0(a),  45°(b),  and  90°(c), 

ka  =  0.  5,  and  different  values  of  p.  * 

Graph  1-8:  |  g+(cp,  cpQ)  | 3  versus  cp  for  cpQ  =  0(a),  45°(b),  and  90°(c), 

ka  =  0.  7,  and  different  values  of  p .  * 

Graph  1-9:  versus  cp*  for  cpQ  =  0(a),  45°(b),  and  90°(c), 

ka  =  0.  9,  and  different  values  of  p.  * 

Graph  1-10:  |  g+(<p>  <pQ)  I  3  versus  cp  for  cpQ  =  0(a),  45°(b),  and  90  (c) 

ka  =  1.1,  and  different  values  of  p.  * 

❖  The  number  next  to  a  curve  gives  the  value  of  p  =  b/a  = 

(minor  axis/major  axis). 
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Graph  1-5: 


Graph  1-6: 
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Graph  1-2  (Con’t):  |g+( 
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Graph  I-6(b):  |  g+(cp,  45°)  |  3  . 
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Graph  I-8(a):  |g+(«p.  0)  |  8  . 
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Graph  1-8 (b) :  |g+(cp,  45°)  | 
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Graph  I-9(b):  |  g+(cp,  45°)  |  2 
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Graph  I-9(c):  |  g+(cp,  90°)  |  » 
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Graph  I-10(b):  J  g+(ep,  45°)  |  ® 
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SET  II:  Graphical  Results  for  Elliptic  Cylinders  and  E  Parallel. 


Graph  II- 1:  -Reg  (cpQf  cpQ)  versus  cpQ  for  ka  =  0.5,  0.7,  0.9, 

and  1.1,  and  different  values  of  p.  * 

Graph  II-Z:  |g  (cpQ,  cp  Q)  |  3  versus  cpQ  for  ka  =  0.5,  0.7,  0.9,  and  1.1 

and  different  values  of  p.  * 

Graph  II- 3:  I  Vq)  I  2  and  the  phase  of  g_(cpQ>  cpQ)  versus  ka^l.l, 

for  cpQ  =  0°(a),  45°(b),  and  90°(c),  and  different  values  of  p.  * 

Graph  II-4:  |g_(iH-ep0,  cp0)  |  2  versus  cpQ  for  ka  =  0.  5,  0.  7,  0.  9,  and  1.  1 

and  different  values  of  p.  * 

Graph  II- 5 :  |  g_(Tr+cp0,  cpQ)  |  s  and  the  phase  of  cpQ)  versus  ka<l.l, 

for  cpQ  =  0°(a),  45°(b),  and  90°(c),  and  different  values  of  p.  * 


*  t 


Graph  II-6: 

|g_(cp>  <P0)  1 

2  versus  cd  for  cpQ 

= 

0(a), 

45°(b), 

and  90°(c), 

ka  =  0.3, 

and  different  values 

of 

p.* 

Graph  II- 7: 

1  g_(<p-  <p0)  1 

2  versus  cp  for  cpQ 

= 

0(a), 

45°(b), 

and  90°(c), 

ka  =  0.5, 

and  different  values 

of 

D# 

Graph  II- 8 “ 

|g.(«p*  «p0)  1 

2  versus  cp  for  cpQ 

= 

0(a), 

45°(b), 

and  90°(c), 

ka  =  0.  7, 

and  different  values 

of 

P-* 

Graph  II-9: 

|g.(«P>  (PoH 

2  versus  cp  for  cpQ 

= 

0(a), 

45°(b), 

and  90°(c), 

ka  =  0.  9, 

and  different  values 

of 

P-* 

> 

Graph  II- 10 

IgJv*  90)  1 

3  versus  cp  for  cpQ 

= 

0(a), 

45°(b), 

and  90°(c), 

I 

.  i 

} 

ka  =  1.1, 

and  different  values 

of 

P  . 

*The  number  next  to  a  curve  gives  the  value  of  p  =  b/a  = 
(minor  axis /major  axis). 
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Graph  II- 1:  -Reg  (epQ,  ?q)* 
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Graph  II-l  (Con't):  -Re  g_(cpQ,  c pQ)  - 
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Graph  II-4:  |g_(-rr+cp 
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ka=0. 5 


Graph  II- 7(b):  '  jg_(cp,  45°) 


l» 


ka=0.  5 
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Graph  II-7(c) 


g-(cp.  90°)  | a  . 
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Graph  II-9(b);  |  gjcp,  45°)  |  *  . 
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Graph  II-9(c) 


DIRECTION  OF 
I  NC I DENCE 


ka=l.  1 


Graph  II-l 0(b):  | g_(cp,  45°)  | a  . 
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Graph  II-l  0(c):  |  g_(cp,  90°)  |  a  . 
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SET  III:  Graphical  Results  for  Semi-Elliptic  Protuberances  and  H  Parallel. 

Graph  III- 1:  -Re  f+(cpQl  TT-cp0)  versus  cpQ  for  ka  =  0.5,  0.7,  0.9, 

and  1.1,  and  different  values  of  6.  * 

Graph  III-2:  |  f+(cp0,  ir-cpo)  |  2  versus  cpQ  for  ka  =  0.  5,  0.  7,  0.  9, 

and  1.1,  and  different  values  of  6.  * 

Graph  HI-3:  I  f+(<P0>  ^-<P0)  |  3  and  the  phase  of  f+(cpQ,  Tr-cpQ)  versus  ka<l.l, 

for  c pQ  =  0°(a)  and  45°(b),  and  different  values  of  6.  * 

Graph  III-4:  |*+("<P0»ir“90)  1 3  versus  cpQ  for  ka  =  0.5,  0.7,  0.9, 

and  1.1,  and  different  values  of  6.  * 

Graph  III- 5:  |  f_|_(-cp0i  Tr-cpo)  I  2  and  the  phase  of  Tr~cPc>)  versus  ka<l.l, 

for  cp0  =  45°,  and  different  values  of  6.  * 

Graph  III -6:  | f+(cp,  "rr-cpo)  | 2  versus  <p  for  cpQ  =  0,  45°,  and  80°, 

ka  =  0.  3,  and  different  values  of  6.  * 

Graph  III- 7:  |  f+(cp,  TT-cpQ)  |  a  versus  cp  for  cpQ  =  0,  45°,  and  80°, 

ka  =  0.  7,  and  different  values  of  6.* 

Graph  III-8:  |  £+(cp,  Tr-cpQ)  1 3  versus  cp  for  =  0,  45°,  and  90°, 

ka  =  1.1,  and  different  values  of  6 .  * 


*The  number  next  to  a  curve  gives  the  value  of  6  =  7|/ §  = 

(axis  ||  ground  plane) /(axis  ground  plane);  6  =  0  corresponds  to 
perpendicular  strips,  6  =  1  to  semi-circles,  and  6  =  oo  to  flat  strips 
(see  Fig.  3). 
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Graph  III-l  (Con't):  -Re  f  ■  (cp 
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Graph  III- 
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Graph  III- 2  (Con't):  |f  (ep, 
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Graph  III-3(a) 
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Graph  III-4:  |f 
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Graph  III-4  (Con't): 
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9  2 
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Graph  III-4  (Con't):  |f+(-cp0<  ir-cpQ)  I 
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(p  —  DEGREES 

Graph  III- 7  :  ka  =  0.  7 
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Graph  III-8  :  ka  =  1.  1 
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Graph  III-8  (con't)  :  ka  =  1. 1 
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<p  —  DEGREES 

Graph  III-8  (con't)  :  ka  = 

102 


SET  IV:  Graphical  Results  for  Semi-Elliptic  Protuberances  and  g  Parallel. 


Graph  IV- 1:  -Re  f_(cp0,  Tt-cp  )  versus  cp  for  ka  =  0.5,  0.7,  0.9, 


Yo'  tc 

and  1.1,  and  different  values  of  6.  * 


Graph  IV-2:  I | 3  versus  cpQ  for  ka  -  0.5,  0.7,  0.9, 

and  1.1,  and  different  values  of  6.* 


Graph  IV-3:  I  2  anc*  the  phase  of  (cpQ,  tt* —cp0)  versus  ka<l.l, 

for  cp  =  0°(a)  and  45°(b),  and  different  values  of  6.  * 


Graph  IV-4:  |f  (-cpQ,  ir-cpo)  | s  versus  cp  for  ka  =  0.5,  0.7,  0.9, 

and  1.1,  and  different  values  of  6.  * 


Graph  IV-5:  |  f_(-<pQ,  Tr-cpQ)  |  3  and  the  phase  of  f__(-cpo>  TT"*^P0)  versus  ka^l.l, 

for  cp  =45°  and  different  values  of  6.* 

'  o 

Graph  IV-6:  |  £_  (cp,  it  -  cpQ)  |  2  versus  cp  for  cpQ  =  0,  45°,  and  80°, 

ka  =  0.  3,  and  different  values  of  6.  ❖ 


Graph  IV- 7:  |  f _ (cp  ,  7r-cp0)  |  2  versus  cp  for  cpQ  =  0,  45°,  and  80°, 

ka  =  0.  7,  and  different  values  of  6.  * 


Graph  IV-8:  |  f_(cp,  Tr”<^P0)  1 2  versus  cp  for  cpQ  =  0,  45°,  and  90°, 

ka  =  1.1,  and  different  values  of  6.  * 


❖  The  number  next  to  a  curve  gives  the  value  of  6  =  T|/ 5  = 

(axis  ||  ground  plane)/ (axis  ground  plane);  6  =  0  corresponds  to 
perpendicular  strips,  6  =  1  to  semi-circles,  and  6  =  00  to  flat  strips 

(see  Fig.  3). 
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<Po  =  90' 


ka  =  0.  9 


Graph  IV- 1  (Con't):  -Re  f  (cp  ,  ir-cp  ). 

TOO 
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(O  =  90° 


cpn  =  90° 


Graph  IV-4  (Con't):  |f  ( -cpQ , tt - cp^)  | 3  . 


1  1  1 
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0 . 04- 


Graph  IV-6  (con't)  :  ka  =  0.  3 
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